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Algebraic reduction for space-time codes 
based on quaternion algebras 

L. Luzzi G. Rekaya-Ben Othman* J.-C. Belfiore 



Abstract 



CD ' In this paper we introduce a new right preprocessing method for the decoding of 2 x 2 algebraic 

o . 

CO , of the new reduction is to absorb part of the channel into the code, by approximating the channel matrix 



STBCs, called algebraic reduction, which exploits the multiplicative structure of the code. The principle 



with an element of the maximal order of the algebra. 

We prove that algebraic reduction attains the receive diversity when followed by a simple ZF detection. 
q | Simulation results for the Golden Code show that using MMSE-GDFE left preprocessing, algebraic 

reduction with simple ZF detection has a loss of only 3 dB with respect to ML decoding. 

| Index Terms: Algebraic reduction, right preprocessing, Golden Code 

^ . EDICS category: MSP-DECD 

CO 

ON ( 

I. Introduction 

OO 
O 



Space-time coding for multiple antenna systems is an efficient device to compensate the effects of 
fading in wireless channels through diversity techniques, and allows for increased data rates. 
A new generation of space-time code designs for MIMO channels, based on suitable subsets of division 
algebras, has been recently developed 11711 . The algebraic constructions guarantee that these codes are 
full-rank, full-rate and information-lossless, and have the non-vanishing determinant property. 
Up to now, the decoding of algebraic space-time codes has been performed using their lattice point 
representation. In particular, maximum likelihood decoders such as the Sphere Decoder or the Schnorr- 
Euchner algorithm are currently employed. However, the complexity of these decoders is prohibitive for 
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practical implementation, especially for lattices of high dimension, arising from MIMO systems with a 
large number of transmit and receive antennas. 

On the other side, suboptimal decoders like ZF, DFE, MMSE have low complexity but their performance 
is poor; in particular they don't preserve the diversity order of the system. 

The use of preprocessing before decoding improves the performance of suboptimal decoders, and reduces 



considerably the complexity of ML decoders 11211 . Two types of preprocessing are possible: 

- Left preprocessing (MMSE-GDFE) to obtain a better conditioned channel matrix; 

- Right preprocessing (lattice reduction) in order to have a quasi-orthogonal lattice. The most widely 
used lattice reduction is the LLL reduction. 

We are interested here in the right preprocessing stage; we propose a new reduction method for 2 x 2 
space-time codes based on quaternion algebras which directly exploits the multiplicative structure of the 
space-time code in addition to the lattice structure. Up to now, algebraic tools have been used exclusively 
for coding but never for decoding. Algebraic reduction consists in absorbing a part of the channel into the 
code. This is done by approximating the channel matrix with a unit of a maximal order of the quaternion 
algebra. 



The algebraic reduction has already been implemented by Rekaya et al. II16H for the fast fading channel, 
in the case of rotated constellations based on algebraic number fields. In this context, the units in the ring 
of integers of the field form an abelian multiplicative group whose generators are described by Dirichlet's 



unit theorem [10]. The reduction algorithm then amounts to decoding in the logarithmic lattice of the 
unit group, which is fixed once and for all. In this case, one can show that the diversity of the channel 
is preserved. 

For quaternion skewfields, which are the object of this paper, the situation is more complicated because 
the unit group is not commutative. However, it is still possible to find a finite presentation of the group, 
that is a finite set of generators and relations. 

Supposing that a presentation is known, we describe an algorithm to find the best approximation of the 
channel matrix as a product of the generators. 

As an example, we consider the Golden Code, and find a set of generators for the unit group of its maximal 
order. Our simulation results for the Golden Code show that using MMSE-GDFE left preprocessing, the 
performance of algebraic reduction with ZF decoding is within 3 dB of the ML. 

The paper is organized as follows: in Section JI] we introduce the system model; in Section HIT] we 
explain the general method of algebraic reduction. In Section JV] we present the search algorithm to 
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approximate the channel matrix with a unit; in Section |V] we prove that our method yields diversity 
order equal to 2 when followed by a simple ZF decoder. We discuss its performance obtained through 
simulations in the case of the Golden Code, and compare algebraic reduction and LLL reduction using 
various decoders (ZF, ZF-DFE), with and without MMSE-GDFE preprocessing. Finally, in Section [Vj 
we describe a method to obtain the generators of the unit group for a maximal order in a quaternion 
algebra. The computations are carried out in detail for the case of the algebra of the Golden Code. 

II. System model and notation 

A. System model 

We consider a quasi-static 2x2 MIMO system employing a space-time block code. The received 
signal is given by 

Y = HX + W, X,H,Y, W G M 2 (C) (1) 

The entries of H are i.i.d. complex Gaussian random variables with zero mean and variance per real 
dimension equal to \, and W is the Gaussian noise with i.i.d. entries of zero mean and variance iVo. 
X is the transmitted codeword. In this paper we are interested in STBCs that are subsets of a principal 
ideal Oa of a maximal order O in a cyclic division algebra A of index 2 over Q(i) (a quaternion algebra). 



We refer to 11711 for the necessary background about space-time codes from cyclic division algebras, and 



to |8|] for a discussion of codes based on maximal orders. 

Example (The Golden Code). The Golden Code falls into this category (see ill] and ||8|]). It is based on 
the cyclic algebra A = (Q(i, 9)/Q(i),a,i), where 9 = and a : x h-> x is such that a{9) = 9 = 1 — 9 
and a leaves the elements of Q(i) fixed. 
It has been shown in 18|] that 

0= \ I Xl X2 J , x 1) x 2 eZ[i i 9] \ (2) 

[ \ 1X2 X X ) J 

is a maximal order of A. O can be written as O = Z[i, 9] © Z[i, 9]j, where 

'-(::) » 

Up to a scaling constant, the Golden Code is a subset of the two-sided ideal Oa = aO, with a = 1 + i9 



llllll . Every codeword of Q has the form 

X = 
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with x\ = s\ + S2O, X2 = S3 + S4#. The symbols si, S2, S3, S4 belong to a QAM constellation. 
B. Notation 

In the following paragraphs we will often pass from the 2x2 matrix notation for the transmitted and 
received signals to their lattice point representation as complex vectors of length 4. To avoid confusion, 
4x4 matrices and vectors of length 4 are written in boldface (using capital letters and small letters 
respectively), while 2x2 matrices are not in bold. 

Notation (Vectorization of matrices). Let <fi be the function M 2 (C) — > C 4 that vectorizes matrices: 

a c\ . 

\^(a,b,c,d) t (4) 

b dj 

The left multiplication function Ai : M 2 (C) — > M2(C) that maps B to AB induces a linear mapping 
A l = (j) o Ai o ^T x : C 4 C 4 . That is, 

<j){AB) = Ai(p(B) VA, B G M 2 (C) 

A; is the block diagonal matrix 

(a o\ 

A; = (5) 

Notation (Lattice point representation). Let {wi,W2,ws, w^} be a basis of aO as a Z[i]-module. Every 
codeword X can be written as 

4 

X = J ^2s i w i , s = (si,S2,s 3 ,s 4 )* e Z[i] 4 
i=i 

Let $ be the matrix whose columns are 

<j)(wi),(j)(w2),(p(w 3 ),(p(w4) (6) 

Then the lattice point corresponding to X is 

4 

x = <p(X) = ^ s i<P( w i) = * s 
i=i 

We denote by A the Z[i]-lattice with generator matrix 

A complex matrix T is called unimodular if the elements of T belong to and det(T) G {1, — 1, i, —i}. 
Recall that two generator matrices $ and span the same -lattice if <J>' = <I>T with T unimodular. 
The following remark explains the relation between the units of the maximal order O of the code algebra 
and unimodular transformations of the code lattice. This property is fundamental for algebraic reduction. 
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Remark 1 (Units and unimodular transformations). Suppose that U G O* is an invertible element: 
then {Uwi, Uw2, Uws, Uw^} is still a basis of aO seen as a Z[z] -lattice. The codeword X can also be 
expressed in the new basis: 

4 

X = ^2 s'i(Uwi), s' = (si, 4) s 3! 4)* 6 Z[i] 4 
i=i 

The vectorized signal is 

4 4 4 

i=l i=l i=l 

Now consider the change of coordinates matrix Try = $~ 1 U^ G M^(C) from the basis {^>(n?i)}i=i,...,4 
to {4>(Uwi)}i=i t ,,,^. We have det(Tfy) = det(U') = det(C/) 2 = ±1, see equation ([5]). Moreover, we 
have seen that Vs G Z[i] 4 , s' = Tf/s G Z[i] 4 . Then Try is unimodular, and the lattice generated by <&Try 
is still A. 

III. Algebraic reduction 

In this section we introduce the principle of algebraic reduction. First of all, we consider a normalization 
of the received signal. In the system model £[]), the channel matrix H has nonzero determinant with 
probability 1, and so it can be rewritten as 



H = y/dsb(H)Hu #i G SL 2 {C) 

Therefore the system is equivalent to 

Y\ = ; Y = HiX + W x 

Algebraic reduction consists in approximating the normalized channel matrix Hi with a unit U of norm 
1 of the maximal order O of the algebra of the considered STBC, that is an element U of O such that 

det(C7) = 1. 

A. Perfect approximation 

In order to simplify the exposition, we first consider the ideal case where we have a perfect approxi- 
mation: H\ = U. Of course this is extremely unlikely in practice; the general case will be described in 
the next paragraph. 
The received signal can be written: 

Y 1 = UX + Wi (7) 
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and UX is still a codeword. In fact, since U is invertible, 

{UX | X G Oa} = Oa 

Applying <f> to both sides of equation ©, we find that the equivalent system in vectorized form is 

yi = U ; *s + wi 

where <& is the matrix defined in ©, s G Z[i] 4 , ji = 4>{Yi), wi = 4>{W{). 
We have seen in Remark Q] that since U is a unit, 

with Tjy unimodular. So 

yi = *Tf/s + wi = *si + wi, si G Z[i] 4 
In order to decode, we can simply consider ZF detection: 



§1 = [* _1 yi] 



81 + 



1 



:<!> W 



where [ ] denotes the rounding of each vector component to the nearest (Gaussian) integer. 

If $ is unitary, as in the case of the Golden Code, algebraic reduction followed by ZF detection gives 

optimal (ML) performance. 



B. General case 

In the general case, the approximation is not perfect with probability 1 and we must take into account 
the approximation error E. We write H\ = EU, and the vectorized received signal is 

yi = E;U;$S + wi = Ez*TVs + Wi = E;*Si + Wi 

The estimated signal after ZF detection is 



si = [S^Ef Vi] 



si + 



=* -1 Er 1 w 



[si + nj 



(8) 



Finally, one can recover an estimate of the initial signal s = T f/ 1 §i. 



Thus, the system is equivalent to a non-fading system where the noise n is no longer white Gaussian. 
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C. Choice of U for the ZF decoder 

We suppose here for simplicity that the generator matrix $ is unitary, but a similar criterion can be 
established in a more general case. We have seen that ideally the error term E should be unitary in order 
to have optimality for the ZF decoder, so we should choose the unit U in such a way that E = H\U~ l 
is quasi-orthogonal. We require that the Frobenius norm \\E\\ F should be minimized^: 

U = argmm\\UH{ l \\ 2 (9) 

ueo, 

det(f/)=l 

This criterion corresponds to minimizing the trace of the co variance matrix of the new noise n in ([8]): 
Cov(n) = Gov ( 7 =^=*-E-w) = J^y^ Cov(w) (E-)* = 



N 



\det(H)\ 



S^Ef 1 (E^f 



and 



IV. The approximation algorithm 

In this section we describe an algorithm to find the nearest unit U to the normalized channel matrix 
Hi with respect to the criterion (©. To do this we need to understand the structure of the group of units 
of the maximal order O. 

Notation. We denote elements of ST^C) with capital letters (for example H\, U) when considering their 
matrix representation, and with small letters (for example h±,u) when we want to stress that they are 
group elements. 

Remark 2 (Units of norm 1). The set 

O 1 = {u G O* | det(u) = 1} 

is a subgroup of O. 

In fact, if u is a unit of the Z[i]-order O, then ^A/Q(i)( u ) = det(it) is a unit in that is, det(u) G 
{1, — 4'}. O 1 is the kernel of the reduced norm mapping N = ^a/QU) '■ ®* ~~ ^ {!> — 1> *i — which 
is a group homomorphism, thus it is a subgroup of O. 

'Remark that since det(E) = 1, \\E\\ 2 F = 
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Table I 

Generators of O 1 . 



id 

U\ = | | =W, u, - = | | = lb 

id 

u 2 = | I = i + (1 + i)j u 2 = I J = i — (1 + i)j 

, -i + 1 i 

( e , / 

U3 = 



= 6>-(l + i)j uj 1 = 

= (l+i) + (l+*0)j u^ 1 = I i! " ^ I =(l + i) + (l + iS)j 

\- i(l + i0) l+i J 

I l+i -l-i9\ 
(1 + i) + (1 + i8)j ue 1 ^ _ = (l + i) - (l + M)j 



lie, 



i6 


\ 

= i9, 


id 


I 


i 




i - 1 

e 


i J 
l+i\ 


i - 1 

e 


9 J 

"'<"') 


-i + l 




l + i 


l + i8 



J(l + i6) l + i 
l + i 1 + it 
i(l + i9) l + i 



ur=\ 1 * 9 + * I = (1 - t) + (9 + i)j u? 1 = I * " " " I = (1 - i) - (9 + i)j 

\i(8 + i) l-i 

l T l 6 + 1 I =(l-i) + (9 + i)j Wg 1 = I ' " " | =(l-i)-(9 + i)j 

i(6 + i) l-i I \-i(9 + i) l-i 



9 




-i + l 


6 J 


9 1 


+ i\ _ 


- 1 


6 ) 


l + i 


-1 - it 


-i(l + it 


) l + i 


l+i 


-1-it 


+ lb 


) l + i 


1 - i 


-9-i\ 


-i(9 + i) 


l-i J 


1 - i 


-6-i\ 



Example (The Golden Code). In the case of the Golden Code, N is surjective since N(l) = 1, N(0) = 
99 = -1, N(j) = -j 2 = -i, N(j9) = i. So {1, -1, i, -i} ^ 0*/O l , and O 1 is a normal subgroup of 
index 4 of 0*. In order to obtain a set of generators, it is then sufficient to study the structure of O . 
Its cosets can be obtained by multiplying for one of the coset leaders {l,6,j,9j}. 

Our problem is then reduced to studying the subgroup O l . In particular, we need to find a presentation 
of this group: a set of generators S and a set of relations R among these generators. In fact, one can 
show that O 1 is finitely presentable, that is it admits a presentation with 5 and R finite (9D. 

Example (Generators and relations in the case of the Golden Code). The group O 1 is generated by 
8 units, that are displayed in Table HI The corresponding relations are shown in Table |II1 



The method for finding a presentation is based on the Swan algorithm 11811 . As it is not well known 
and is rather complex, we have chosen to expose it in detail in Section ED 
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Table II 

Fundamental relations among the generators of O 1 . 



u| = -l 

ui = -i 

(lllUlf = 1 

(U2U- 1 ) 3 = 1 

U6M3U7 = —1 

M6W7W4 1 = —1 

U8U3U5 = — 1 

Mj 1 M8M5 = —1 

WlttJ MlMj 1 = 1 

U^ 1 U2U^ 1 U^ 1 UsU2UsUi = 1 

U6M^" 1 W6UlW7~ 1 M^" 1 W^" 1 U 1 " 1 = 1 



A. Action of the group on the hyperbolic space H 3 

The search algorithm is based on the action of the group on a suitable space. We use the fact that O l 
is a subgroup of the special linear g roup ST^C), and consider the action of SX2(C) on the hyperbolic 
3-space H 3 (see for example |[5|] or H 1 3M for a reference). 

We refer to the upper half-space model of H 3 : 

M 3 = {(z, r) | z G C, r G M, r > 0} (1 1) 

H 3 can also be seen as a subset of the Hamilton quaternions H: a point P can be written as (z, r) = 
z + r j = x + iy + rj, where {1, i, j, k} is the standard basis of TL. We endow H 3 with the hyperbolic 
distance p such that if P = z + rj, P' = z' + r'j, 

coshp(P,P0 = l+ d( ^ /)2 , 

where d(P,P') 2 = \z — z'\ 2 + (r — r') 2 is the squared Euclidean distance. The corresponding surface 
and volume forms on H 3 are ( 01311 . pp. 48-49) 

ds= dx2 + d ^ + dr \ (12) 

do = (13) 
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The geodesies with respect to this metric are the (Euclidean) half-circles perpendicular to the plane 
{r = 0} and with center on this plane, and the half-lines perpendicular to {r = 0}. Given a matrix 

its action on a point P = (z, r) is defined as follows: 

_ (az+b)(cz+d)+acr 2 

\ cz+ d\<+\c\^ ' (U) 

# T 

r - \cz+d\ 2 + \c[ J r 2 

(Here we denote by z the complex conjugate of z). 

The action of g and —g is the same, so there is an induced action of PSL2(C) = 5^2 (C) /{I, — 1}. 
PSL2(C) can be identified with the group 



som + (BI 3 ) of orientation-preserving isometries of H 3 with 
respect to the metric defined previously ( 111311 . p. 48). 

All the information we will gain about the group O 1 will thus be modulo the equivalence relation g ~ — g; 
we denote by PO l its quotient with respect to this relation. 
Consider the action of PSL2(C) on the special point 

J = (0,l)=j (15) 

which has the following nice property (|5J, Proposition 1.7): 

Vg£SL 2 (C), \\g\\ 2 F = 2 cosh p(J,g(J)) (16) 

Remark 3. If g £ U(2) is unitary, then g leaves every point of H 3 fixed ([5], Proposition 1.1). Then by 
considering for example the mapping PSL2(C) — > H 3 that sends g to g(J), one can identify H 3 with 
the quotient space PSL 2 (C)/U(2). 

B. The algorithm 

We assume here the following fundamental properties, which will be proven in Section [VTJ 

1) {u(J) | u £ O 1 } is a discrete set in H 3 . 

2) Given a unit bGO 1 , the set 

V u = {P e M 3 | p{P, u(J)) < p(P, u'{J)) Vu' / u} 
is a compact hyperbolic polyhedron with finite volume and finitely many faces. 
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3) Two distinct polyhedra V U ,V U ' can intersect at most in one face; all the polyhedra are isometric, 
and they cover the whole space H 3 , forming a tiling. Moreover, if V = Vt is the polyhedron 
containing J, 

V u = u(V) 

4) The polyhedra adjacent to V are given by 

U1 (V), u r (V),u^(V), u-\V) (17) 

where {u±, . . . , u r } is a minimal set of generators for O l . 
As anticipated in Section [TTll given the normalized channel matrix h\ G SL2 (C) we want to find 



u = argmin \\uh 1 1 || p (18) 

u&O 1 



But we know from equation (fT6l ) that 



Hn/i- 1 )^ = 2cosh(/j(J,n/i^ 1 (J))) = 2cosli(p(ir 1 (J) ) h^(J))), 

since u is an isometry. So the condition (fT8l) is equivalent to 

n = argmin J), /i^ 1 (J)) 

The point h^ l (J) is contained in the image n(P) = 7^ of the polyhedron V for some u G O 1 . It follows 
from the definition of Vu that h^ l (J) is closer to u(J) than to any other u(J), u G O 1 . Since all the 
polyhedra are isometric, 

max 

where i? max is the radius of the smallest (hyperbolic) sphere containing V. Therefore we have the 
following property: 

\uK{ x f F <C (19) 



We now go back to the problem of finding a unit u such that h\ (J) G u(V), given a normalized 
channel matrix hi G SI^C). Let u\, . . . , u r be the generators of O 1 in (fT71) and u r+ i = u^ 1 , . . . , U2 r = 
u~ l their inverses. The neighboring polyhedra of V are all of the form Ui(V), i = 1, . . . , 2r. 
The idea is to begin the search from V and the neighboring polyhedra, corresponding to the generators 
of the group and their inverses, and choose the Ui such that Ui(J) is the closest to h^ 1 (J). Since is 
an isometry of H 3 , at the next step we can apply u^ 1 and start again the search of the Ui> that gives the 
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Figure 1. A step of the algorithm. The polyhedra are represented as two-dimensional polygons for simplicity. 

closest point to Ui~ 1 /iJ~ 1 (J). With this strategy we only need to update a single point and perform 2r 
comparisons at each step of the search. 

The algorithm is illustrated in Figure [TJ 
Suppose that the matrix form of the has been stored in memory at the beginning of the program, 
together with the images u\(J),... ,U2 r {J) of J, for example using the coordinates in the upper half- 
space model (fTTb . Let 

Ui(J) = (xi,yi,n), i = l,...,2r 

INPUT: hi G 5L 2 (C). 

Initialization: let h = h\,u = 1,zq = 0. 

REPEAT 

1) Compute h~ l (J) = (x,y,r). 

2) Compute the distances 

dk = 2 coshp^-V), = 1 + (X ~ X * )2 + ( " )2 + (r ~ r * )2 ' i = 1, • • • , 2r, 

2rr; 

do = 2coshp(/i~ 1 (J), J) 

3) Lef io = argmirij 6 ro j x i ... i 2r-} (If several indices i attain the minimum, choose the smallest.) 

4) Update u <— ?Ziij , /i <— /mj . 
UNTIL t = 0. 

OUTPUT: u = u^ 1 is the chosen unit. 
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Remark 4 (Advantage of the algebraic reduction in the case of slow fading channels). If the 

channel varies slowly from one time block to the next, it is reasonable to expect that the polyhedron 
u(V) containing (J) at the time t will be the same, or will be adjacent, to the polyhedron chosen at 
the time t — 1. Thus, this method requires only a slight adjustment of the previous search at each step. 
On the contrary, the LLL reduction method requires a full lattice reduction at each time block. 

V. Performance of the algebraic reduction 

A. Diversity 



It has recently been proved 11911 that MIMO decoding based on LLL reduction followed by zero-forcing 
achieves the receive diversity. The following Proposition shows that algebraic reduction is equivalent to 
LLL reduction in terms of diversity for the case of 2 transmit and 2 receive antennas: 

Proposition 1. The diversity order of the algebraic reduction method with ZF detection is 2. 

Proof: We suppose that the symbols Si, i = 1,...,4 belong to an M-QAM constellation, with 
M = 2 2m . Let £ av be the average energy per symbol, and 7 = ^ the SNR. 

For a fixed realization of the channel matrix H, equation ([8]) is equivalent to an additive channel without 
fading where the noise n is no longer white. 

We can compute the error probability using ZF detection conditioned to a certain value of H, and then 
average over the distribution of H: 

P e ( 7 ) = J Pe{l\ H)dH (20) 

With symbol by symbol ZF detection, P e (7) is bounded by the error probability for each symbol: 

4 



i=l 



Using the classical expression of P e in a Gaussian channel, for square QAM constellations ( 11511 . §5.2.9), 
we obtain ( 

3£ a 



P((si)i + < 4erfc I y g2 ( M _ ^ J ^ 4e ^ 

where er 2 is the variance for complex dimension of the noise component rii . 

We have seen in ( fTOl that the trace of the covariance matrix of the new noise n is bounded by 
| dc ^jj-) | ||* ll-^r 1 !!^' recamn § tnat the Frobenius norm is submultiplicative. Thus 

o?<Cov(n)< CN ° 



|det(F)| 
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2 i 1 1 2 

because ^ = ||f7-fff \\ F < Co, see equation ( fl9l ). 

Indeed if is unitary, as in the case of the Golden Code, ||$ _1 E i _1 ||^, 

Finally, 



e; 



-in 2 



< Co. 



P e ( 7 | H) < 16e~ 



)|dct(H)H 



16e" c|det(if)l7 



In order to compute the error probability in equation d20b . we need the distribution of |det(iT)|. It is 
known 17i|4|] that if H is gaussian with i.i.d. M(0, 1) entries (variance per real dimension |), the random 
variable 4 |det(i?)| 2 , corresponding to the determinant of the Wishart matrix 2HH H , is distributed as the 
product of two independent chi square random variables with 2 and 4 degrees of freedom respectively. 
Consider two random variables X ~ x 2 (2),^ ~ X 2 (4): their joint probability distribution function is 

Px,y(x, y) = gV e ~*~ * x,y > 
Then the cumulative distribution function of Z = 2 |det(-ff)| = vXY is 



F z (z) = P{VXY <z}= / / Vx,y{x, y)dxdy 

J J ^/xy<z 

From the invertible change of variables u = y, v = yjxy with Jacobian 

J = 

we obtain 



u 2 u 



1 







2v 
u 



F Z {z) 
Pz(z) 



/ v 2 



u 

oo 



PX,Y 
JO 

dF z (z) z 



u I \J\ dudv 



• du 



e 



o 



zdu] dv 



e 2., 



dz 4 J 2 
where J^i is the modified Bessel function of the second kind. Finally, 



Peiri) < E 

= 16 





= 16 [ 




Jo 



1 



+ 



— T 



where ^ is the Digamma function. The series in the last expression being uniformly bounded for large 



7, the leading term is of the order of -ij. 



□ 
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B. Some remarks about complexity 

The length of the algorithm described in Section IIV-BI is related to the initial distance 



2coshp(/i^ 1 (J), J) = \\hl 



in 2 



\h 



i\\f 



H 



VdetH 



\H\ 



\det(H)\ 



In order to have more information about the distribution of this distance, one has to find the distribution 



of the random variable 



\det(H)\- 



From IJ], we learn that H is unitarily similar to 



H 



1 (X 

Y Z 



where X 2 ~ x 2 (4), Y 2 , Z 2 ~ x (2) and X, Y, Z are independent. Therefore 

2 



I -Hi 



H 



X 2 + Y 2 + Z< 



|det(fl)| 



det(F) 



We want to find the distribution of the random variable T 

X, Y, Z: 



XZ 



X2+ xz +Z2 knowing the distributions of 



x 



px{x) 



-e 2 



pz(z) = ze 



Their joint probability distribution is 



Px,Y,z(x,y,z) = -x yze 



.l(x 2 +y 2 +z 2 ) 



and the distribution of T is given by 

PT(t) = ^F T (t) = | 



px,Y,z{x,y,z)dxdydz 



d_ f 00 
dt L 



\/tXZ X 2 "z 2 ~,3„ jr , 2 i 2 i 2 

X yZ _x +y +z 

e 2 dydzdx 

o 2 



With the change of variables % = cosh(u), w = y 2 this integral becomes 

^ / poo />a;(cosh(ii)+sinh(n)) p2cosh(u)xz—x 2 —z 2 ^,3 
<9 W WO A(cosh(ii)-sinh(«)) JO 



X"Z 



dwdzdx 



\ du 



at 



-tanh ( U )_ = _ jn _ 



The following example shows that the distance 1 (J), J) is mostly concentrated near the origin: 

Example. In Section IV"fl we will see that the minimum of the distances between J and the vertices of V 
for the Golden Code is i? m i n = arccosh( 1.9069 •••) = !, 2614 • • • . From the distribution of T, we find 
that the probability that p(J, h^ 1 (J)) > i? m i n is approximately 0.038: in most cases /i^~ 1 (J) is already 
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Figure 2. Performance of algebraic reduction followed by ZF or ZF-DFE decoders using 4-QAM constellations. 



contained in V or in one of the neighboring polyhedra, and the algorithm stops after one step ! Moreover, 
the probability that p(J, J)) > 5R m i n is of the order of 10~ 10 , which is negligible since for practical 
values of the SNR the probabilities of error for ML detection are typically of the order of 10~ 6 at best. 

C. Simulation results 

Figure |2] shows the performance of algebraic reduction followed by ZF and ZF-DFE decoding compared 
with ML decoding using 4-QAM constellations. One can verify that the slope of the probability of error in 
the case of algebraic reduction with ZF detection (without preprocessing) is very close to —2, confirming 
the result of Proposition Q] concerning the diversity order. 



One can add MMSE-GDFE left preprocessing to solve the shaping problem for finite constellations II12H 
in order to improve this performance. With MMSE-GDFE preprocessing, algebraic reduction is within 
4.2 dB and 3.2 dB from the ML using ZF and ZF-DFE decoding, at the FER of 10 -4 . 
In the 16-QAM case, the loss is of 3.4 dB and 2.6 dB respectively for ZF and ZF-DFE decoding at the 
FER of 1CT 3 (Figure In the same figure we compare algebraic reduction to LLL reduction using 
MMSE-GDFE preprocessing. The two performances are very close; with ZF-DFE decoding, algebraic 
reduction has a slight loss (0.3 dB). On the contrary, with ZF decoding, algebraic reduction is slightly 
better (0.4 dB gain), showing that the criterion (O is indeed appropriate for this decoder. 

Numerical simulations also evidence that the average complexity of algebraic reduction is low. In 
Section IIV-BI we have seen that each step of the unit search algorithm requires only a few operations. 
Table [III] shows the actual distribution of the number of steps in the unit search algorithm. The data 
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Figure 3. Comparison of algebraic reduction and LLL reduction using MMSE-GDFE preprocessing combined with ZF or 
ZF-DFE decoding with 16-QAM constellations. 



Table III 

Number of steps of the search algorithm. 



average number of steps 


distribution of the number of steps 




1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


> 11 


1.923062 


38.2% 


39.4% 


16.0% 


4.8% 


1.2% 


0.2% 


3.7 ■ 10~ 2 % 


7.6 • 10~ 3 % 


2.6 • 10~ 3 % 


10~ 4 % 






refers to a computer simulation for the Golden Code using a ZF decoder, for 16-QAM constellations, for 
the transmission of 10 6 codewords. (Clearly this distribution does not depend on the SNR.) The average 
length of the algorithm is less than 2. 

VI. Finding the generators 

In this Section, we describe a method to find a presentation for the group O 1 of units of norm 1, and 
the corresponding polyhedron V. The computations are carried out in detail for the Golden Code. 

A. Kleinian groups and Dirichlet polyhedra 

We introduce some terminology that will be useful later: 

Definition 1 (Kleinian groups). Let Y be a subgroup of the projective special linear group PSZ^C) 
acting on the hyperbolic space H 3 . 

- If T is discrete, that is if the subspace topology on V is the discrete topology, T is called a Kleinian 
group. Remark that then T is countable. 
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- The orbit of a point xq G H 3 is the set {g(xo) \ g E F}. 

- A fundamental set for the action of T is a subset of H 3 containing exactly one point for every orbit. 

- A fundamental domain for T is a closed subset D of H 3 such that 

a) [] geT g(D)=n\ 

b) If g G r \ {1}, the interior of D is disjoint from the interior of g{D). 

c) The boundary of D has measure 0. 

- T is called cocompact if it admits a compact fundamental domain; we say that F has finite covolume 
if it admits a fundamental domain with finite volume. 

- If T has finite covolume, and D\ and D2 are fundamental domains for F, then Vol(-Di) = Vol(Z>2) < 
00 ([13], Lemma 1.2.9). 

In the case of a Kleinian group F, one can obtain a fundamental domain that is a hyperbolic polyhedron 



1131 



2D- This polyhedron can be obtained as an intersection of hyperbolic half-spaces. 
For any pair of distinct points Q,Q' E H 3 , the set of points equidistant to Q and Q' with respect to p 
is a hyperbolic plane, called the bisector between Q and Q', which divides H 3 into two open convex 
half-spaces, one containing Q and the other containing Q'. Given g G F, let 

D g (Q) = {P G M 3 I p(Q, P) < p(g(Q),P)} (21) 

the closed half-space of the points that are closer to Q than to g(Q). If Q is not fixed by any nontrivial 
element of F, the Dirichlet fundamental polyhedron of F with center Q is defined as the intersection of 
all the bisectors corresponding to nontrivial elements: 

Vv = f| D g (Q) (22) 

The definition (1221 cannot be used directly to compute the polyhedron, since we ought to intersect an 
infinite number of bisectors. Let B(Q,R) denote the closed ball with center Q and radius R, and let 

Dr{Q) = f]{D g (Q) \g^l, g(Q) G B(Q, R)} (23) 

If Vr is compact, it has finite diameter, so there exists R > such that Vv = Dr(Q)- 

B. Poincare's theorem 

From the Dirichlet polyhedron of a Kleinian group one can obtain a complete description of the latter, 
including generators and relations. In fact, a famous theorem due to Poincare establishes a correspondence 
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between a set of generators of the group and the isometries which map a face of the polyhedron 
into another face, called side-pairings. The sequences of side -pairings which send an edge into itself 
correspond to a complete set of relations among the generators. 

A complete exposition of Poincare's theorem in the general case can be found in [60. We only need a 
rather weak version of the theorem that we state as follows: 

Theorem 2 (Poincare's polyhedron theorem). Let V be a hyperbolic polyhedron in H 3 with finitely 
many faces. Let T denote the set of faces of V, and suppose that: 

a) [Metric condition] For every pair of disjoint faces ofV, the corresponding geodesic planes have 
no common point at infinity. 

b) [Side-pairings] There exist two maps 1Z : T — > T, U : T — > Isom(H 3 ) such that: 

- VF € T, U 2 {F) = F 

- If 1Z{F) = F', U{F) = uf maps F' onto F, sending distinct vertices into distinct vertices, 
and distinct faces into distinct faces, and maps the interior of V outside of V. Moreover 
u F , = (up)- 1 

1Z is called a side-pairing for V. 

c) [Cycles] For each edge E\ ofV, there is a cycle starting with E\, that is a sequence of the form 
[Ei, . . . , E n+ \], where E{, i = 1, . . . , n + 1 are edges of V, and Vi € {1, . . . , n} there exists a 
generator € U{fF} such that u^(Ei) = Ei+i, and £7 n +i = E\. Moreover, we suppose that 
u = u*™) o • • • o is a rotation through an angle m G Z + , and that its restriction to E\ is 
the identity. 

Consider the group F generated by lA{fF}. Then V is a fundamental domain for the action ofF on H 3 . 



The proof of this theorem is a special case of the proof of Theorem 4.14 in Q6Q . 
C. The structure of O 1 

We now have all the necessary background to find a fundamental domain, and thus a set of generators, 
for PO 1 = O 1 /{I, — 1}. The following theorem shows that PO 1 is a Kleinian group, and describes its 



Dirichlet polyhedron (see |[5|] or Il20ll ): 



Theorem 3. Let A be a quaternion algebra over a number field K such that 

a) K has exactly one pair of complex embeddings 

b) A is ramified at all the real places, that is A ®q K v is a division ring for every real place v of K. 
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Let O be an order of A Then: 

- PO 1 is a Kleinian group. 

- PO 1 has finite covolume and its Dirichlet polyhedron has finitely many faces. 

- PO 1 is cocompact if and only if A is a division ring. 

Remark that conditions (a) and (b) of the theorem are verified since K = Q(i) is an imaginary quadratic 
number field and thus has a pair of complex embeddings and no real embeddings. 

Thus PO 1 admits a compact fundamental polyhedron Vo 1 of the form (|22~1 ). with finitely many faces 
and finite volume. This volume is known a priori and only depends on the choice of the algebra A (see 



I12D, p.336): 



Theorem 4 (Tamagawa Volume Formula). Let A be a quaternion algebra over K such that A <8>j; 
M2(C). Let O be a maximal order of A Then the hyperbolic volume 

VolCPoO = -^Qk{2) \D K \l H(N p - 1) 

p\S c 



In the previous formula, Qk denotes the Dedekind zeta function.] relative to the field K, Dk is the 
discriminant of K, 5q is the discriminant of O, p varies among the primes of Ok, and N p = [Ok '■ pOk], 
where Ok is the ring of integers of K. 

Example (The Golden Code). In the case of the Golden Code algebra, Dq^ = —4, and S(0) = 5Z[i]. 
The only primes that divide the discriminant of the maximal order are (2 + i) and (2 — i), both with 
algebraic norm 5. In conclusion, 

8Com(2)16 32Com(2) 
Vol ?V = W V = = 4, 885149838 • • • (24) 

4tt z ir z 

since Cq(<)(2) = 1.50670301 

Remark 5. We have seen in section ITV-B I that a smaller polyhedron V results in a better average distance 
between u(J) and h~ l {J) and a better approximation. So the algebraic codes such that Vol("P) is small 
are better suited for the method of algebraic reduction. 

2 The Dedekind zeta function is denned as Ck{s) = ^2{[Ok '■ I])", where / varies among the proper ideals of Ok- 

i 
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D. Computing the Dirichlet polyhedron for O 1 

We suppose here that we already have an estimate of the volume of Vo 1 , given by Theorem 4[ For 
this reason our strategy to find the Dirichlet polyhedron differs slightly from that described in yfl. The 
idea is to compute the sequence Dr{Q) defined by (|23T ) for an increasing sequence of values of R, until 
we find R such that D-^{Q) is compact. If the hypotheses of Poincare's Theorem are verified for a set 
of side -pairings belonging to O 1 , D-^(Q) is a Dirichlet polyhedron for some subgroup of O l . To check 
whether this subgroup coincides with O l it is sufficient to estimate of the volume of D-^(Q). 

Following [2], we take as our base point Q the point J defined in ( fT5T ). One needs to check that J is 
not fixed by any nontrivial element of PO 1 . 

Remark 6. As pointed out in [2J], if on the contrary J is fixed by some nontrivial element, one needs 
first to compute a fundamental domain for the stabilizer Fj of J (the subgroup of elements that fix J), 
and then intersect it with p| 5gr \ rj D g (J). 

Because of the property (fT6l ), in order to find Dr(J) we only need to intersect the bisectors corre- 
sponding to elements of O 1 with square Frobenius norm less or equal to 2cosh(i?). Since O can be 
identified with a discrete lattice in C 4 using the map <fi defined in (0]), and the Frobenius norm corresponds 
to the Euclidean norm in C 4 , clearly there is only a finite number of these elements. 

Since cosh is increasing on the positive half-line, in order to find the half-space D g (J) one can solve 

. a b 

the inequality cosh(p(Q, J)) < cosh(p(Q, g( J))). For a general g 

bd + ac 1 



c d , 



9{J) 



|d| 2 + |c| 2 ' M 2 + |c| 2 



and the corresponding half-space has equation 

(C - l)x 2 + (C - l)y 2 + (C - l)r 2 - 2Ax - 2By + - + ^ +1 - 1 > 0, 

where A = $t(bd+ac),B = Ss(bd+ac), C = |d| 2 + |c| 2 . Its boundary is a sphere of center ^^Ay, c?T> ^ 
and square radius ^ f ^jjjp + lj ■ Remark that if we change the sign of the pair a, d or b, c, the radius 
doesn't change, while the center is reflected with respect to the origin. 

If a ball or complementary of a ball (according to the sign of C — 1) in the list (|23l ) is already contained 
in the intersection, we can discard the corresponding element of the group. Since all the spheres have 
center on the plane {r = 0}, in order to determine whether a sphere is contained in another we only 
need to consider their intersections with this plane. 
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Figure 4. The projection of the bisectors on the plane 
{r = 0}. 



Figure 5. The intersection of the spheres (the picture shows 
only half of the space for better understanding). 



E. Computing the generators for the Golden Code 

We now apply the method described in the previous Section to the maximal order O of the algebra of 
the Golden Code. One can easily verify that in this case J is not fixed by any nontrivial element of O . 
Considering the elements g £ O 1 such that \\g\\ 2 F < 9 by computer search, we find that V = D-^(J) is 
compact with R = arccosh (|), since it doesn't intersect the plane "at infinity" {r = 0}. Table U lists 
the elements {ui,u^ l },i = 1, ... ,8 of the group that are necessary to obtain V . The equations of the 
corresponding spheres, the bisectors 

S( Ui ) = D Ui (J), Siu- 1 ) = D u7 i(J), i = l,...,8 

(see definition (|2TI )) can be found in Table [TV] 

The vertices of V are the intersections of all the triples of spheres S(ui): 

Vi, V[ = n'm), V" = *"(Vl), V!" = tt'"^), i = 1 . . . 6, 

Here 7r',7r" and it'" denote the reflections with respect to the plane {y = x}, the plane {y = —x} and 
the line {x = 0, y = 0} respectively. 



v 2 



'5^/5 + 9 3^5- 1 1 



16 



16 



'1 + 



1 9 

2' 2 



, - V 33 + 11VE = S(ui) n 5(u 4 ) n S{u 6 ), 



S(ui)nS(u 2 )nS(u 6 ), 
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Table IV 
Bisectors 



unit 


center 


radius 


int/ext 


Ul 


(0,0,0) 


9 


I 




(0,0,0) 


9 


E 


U2 


(1,1,0) 


1 


E 




(-1,-1,0) 


1 


E 


U3 


(-9,-9,0) 


9 


E 


U3 1 


(9,9,0) 


-9 


E 


U4 


(6,9,0) 


9 


E 


U5 
U(j 

ue 1 
u 7 

1 

Us 

"a" 1 


(-9,-9,0) 

( -9y/5+19 -9-&V5 n\ 

y 22 > 22 ' u j 

(9V5-19 9+5^5 n \ 
^ 22 ' 22 ' U J 

( 9^+19 -9+5v^ n\ 
\ 22 ' 22 ' U J 

/ -9^-19 9-Sv^ q\ 
^ 22 ' 22 ' U J 

( -9-5^/5 -9^5+19 n\ 
\ 22 ' 22 ' U J 

( 9+5^5 9^5-19 r\ 
\ 22 ' 22 ' U J 

( -9+5^5 9^5+19 n\ 
\ 22 ' 22 ' u y 

/ 9-5^5 -9^5-19 n\ 
\ 22 ' 22 ' U ^ 


-9 

22 V 1 V u/ 

^(7-V5) 
^(7 + ^5) 
^(7 + ^5) 
^(7-V5) 

^(7 + ^5) 
4(7 + y/S) 


E 
E 
E 
E 
E 
E 
E 
E 
E 




Figure 6. A schematic representation of the polyhedron Figure 7. The projection of the polyhedron V on the plane 

V (in the picture, the edges have been replaced by straight {r = 0}. 

lines). 
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V * = [-W + 2 ' 10" " 2 ' 2 V 10 J = * M * W ^ } 



The faces of "P correspond to portions F(u) of the spheres S(u), with u one of the units in Table J] 
The projection of the faces of V on the plane {r = 0} is shown in Figure IVI-EI 

As explained in Section IVl-DI in order to prove that V is a Dirichlet polyhedron for O , we will first 
show that it is a fundamental domain for some subgroup T of O 1 using Poincare's Theorem. Comparing 
the volume of V with the value (p4t . we will find that T = O 1 . 

The metric condition in Poincare's Theorem can be verified given the equations of the spheres (see also 
Figure HJ). 

Define a side -pairing as follows: U{F{u)) = u, lZ(F(u)) = F(u _1 ) for every u in Table U The action 
of the generators on the faces and vertices is summarized in Table [Vl] and it is not hard to see that it 
satisfies all the conditions in the theorem. In fact every face F(u) = V n u(V) C S(u). Remark that 
an isometry between polygons with the same number of vertices, sending distinct vertices in distinct 
vertices, must be onto. 

In order to check that the cycle condition of Theorem [2] holds, we need to compute the minimal relations 
or "cycles" between the generators, by finding the sequences of edges of V of the form [E\, . . . , E n+ \], 
such that u^{Ei) = Ei+i, and E n+ i = E\. As must contain Ei, there are only two possible 

choices for «W, corresponding to the two faces containing the edge Ei. 

Given such a sequence, ■ ■ ■ u^> is an element of finite order in O 1 , that is (u^ ■ ■ ■ u^) k = 1 for 
some k. (Remark that every cyclic permutation of the sequence [E\, . . . , E n+ \] gives rise to a new cycle.) 
Actually it is necessary to "lift" the relation from PSL2(C) to 5L2(C). We also require our sequences 
to be irreducible, that is u^ +1 ^ / (u®) -1 for all i. 

In this way we obtain a decomposition of the set of edges of V into cycles. The action of the generators 
on the faces is summarized in Table |VH the cycles are described in Table |V] 

A complete set of relations is listed in Table JI] Except for the first four, the products correspond to the 
identity in PSL2(C) (thus, a trivial rotation). By computing the eigenvalues of ^3,^4,^2^1 and u^u^, 
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Table V 

Cycles for the edges of V. 



ynyn, ynyn, 

-l -l 

V3V4 vi" v-i. vi" vi" -^-> Wv^ 
viVs ^ v 4 "'v 3 '" Wi4 ^ vv 3 

yryT ynyn ynyn ^ yryT 
yryl ^yTyl ^yjfyf ^VTyl 

vvl -^-> Wvf ^ v{"V{' ^ vj"vj' ^ vTv( 
vlyl ^ Wv* ^ VJvJ m> vTv? ^ v[\q^ 

-1 -1 -1 -1 

V1V2 — > VI" V^" V 4 V U -^-> V 4 "'V 6 "' ■ ■ ■ 

VJ/ 6 ^ . . . yTnym yy 2 ^ ynTyTfT T^T 



Table VI 

Action of the generators on the vertices of V. 





')) 


= S(ui) 


m(V 5 ) = V{", ui(Vi) = Vi', «i(W) 


= vi, Ul (Vi") 


= Vi, 








ui(V«) = V 2 '", ui(V 6 ') = V 2 ", ui(V 6 ") 


= Vl «i(V 6 "') 


= Vi 


U2(S(U2 


*)) 


= S(u 2 ) 


U2(Vi) = Vf), u 2 (Vl)=Vl'), 


ua(VZ) = ttf', 










U 2 (V 6 "')=V 2 , U 2 (Vl") = V4, 


u 2 (V 2 "') = V 6 




U3(S(u^ 


')) 


= S(u 3 ) 


u 3 (V s ") = Vi', U3 (Vi") = Vi", u 3 (V B ") 


= W, U8(K") 


= Vi" 




L )) 


= S(u 4 ) 


u 4 (V 3 ) = Va, u 4 {Vi) = Vi, u 4 {Vi) 


= W, ^4(V 5 ) = 


Vy 


u 5 (S(u~ 


1 )) 


= S(u 5 ) 


us(Vi) = Vi', u B (VZ) = V 6 ", U6(V B ') 


= VI', u 5 (Vi) = 


-vi' 




L )) 


= S(u 6 ) 


ue(Vi") = Va, u 6 (Vi") = Vi, u 6 (V{") 


= Va, « 6 (V 2 '") 


= V 2 




L )) 


= S(u 7 ) 


ur(V 3 ) = Kg'", u 7 {V 5 ) = VJ", U7(V 6 ) 


= K", U7(V4) = 


= vi" 


U 8 (S(Ug 


')) 


= 5"(u 8 ) 


^s(V 4 ") = Vi, u a (Vi') = V 2 ', u 8 (Vi') 


= Vi, u s (Vi') -- 


--vi' 
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we find that they are indeed conjugated to rotations of an angle ^ around the axis {x = 0, y = 0}!^ 

We have thus shown that V is a Dirichlet polyhedron for some subgroup T of O 1 . But if T were a 
proper subgroup, the volume of V would be a multiple of the volume of the fundamental polyhedron for 
O 1 that we computed in (El), that is it should be at least 2 • 4.88514 • • • = 9.77029 
So the last step of the proof that V is a fundamental polyhedron for O 1 is the following: 

Lemma 5. Vol(P) < 9.77029 

The proof of this fact is rather tedious and is reported in the Appendix. 

Remark 7. From the coordinates of the vertices of V, one finds that the radius of the smallest hyperbolic 
sphere containing V is 

i?, max = arccosh(2.2360 • • • ) = 1-4436 
while the minimum of the distances between J and the vertices of V is 

R mill = arccosh( 1.9069 • • • ) = 1-2614 • • • 

VII. Conclusions 

In this paper we have introduced a right preprocessing method for the decoding of space-time block 
codes based on quaternion algebras, which allows to improve the performance of suboptimal decoders 
and reduces the complexity of ML decoders. 

The new method exploits the algebraic structure of the code, by approximating the channel matrix with 
a unit in the maximal order of the quaternion algebra. Our simulations show that algebraic reduction 
and LLL reduction have similar performance. However in the case of slow fading, unlike LLL reduction, 
algebraic reduction requires only a slight adjustment of the previous approximation at each time block, 
without needing to perform a full reduction. 

In future work we will deal with the generalization of algebraic reduction to higher-dimensional space 
time codes based on cyclic division algebras. 

3 g £ SL2(C) is called elliptic, and is a rotation around a fixed geodesic, if and only if tr(<?) G R and |tr(i?)| < 2, see Q], 
Prop. 1.4. If its eigenvalues are e' 13 , e _8/3 , then the angle of rotation is 2/3. 
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Appendix 

A. Proof of Lemma \5\ 

In order to prove that the volume of V is smaller than the required constant, we can compute the 
volume of the hyperbolic polyhedron Q enclosed by S(ui), the plane {r = — |}, and the spheres 
5(u 2 ),5(u2 1 ),5(n^ 1 ),5(u 4 ),5(u 3 ), 1 S(uJ 1 ),5(u3 1 ). Clearly Q D V. Recalling the definition of the 
hyperbolic volume in ( fl~3l ), the volume of the spherical sector T enclosed by S{u\) and jr = — 1| 



is 

32 J2 



J _ < 9 ^ r l dr = tt (~ - ln(0) + 2# 4 + In = 36.2937- • • 

To this volume we must subtract the volume of the intersection of T with the chosen spheres. 
From the expression for the area of the intersection of two circles of radii Ri and R 2 whose centers 
have distance d 112 ill 

A(Rx,R 2 , d) = R\ arccos f - + ^ R ~ ^ + R\ arccos f - + ^ R ~ ^ + 



+ 2 V(~ d + Ri- R2)Vi d + Ri- R2){d -Ri + R 2 ){d + Ri + R 2 ), 

we obtain the area of the horizontal sections of T n S(u 2 ). Since R\ = \J Q 2 — f 2 , R 2 = vl — f 2 
are the radii of S(u±) n {r = f}, 5 , (n 2 ) H {r = f } respectively, and the distance between the centers is 
we find 



i? 2 , d) = tt(1 - f 2 ) + (f 2 - 1) arccos ( - A — = ] + 



4 y/l-r 



+ (9 2 — r 2 ) arccos 



a/2^+ 
4 ^ 



±£=) - -V-l + 60 2 -d*-8v 2 , 
2 _ f 2 I 2 



which is defined for r < ^ 9 ^ 3v ^ . In conclusion, 



V9 + 3V/5 



Voi(T n 5(u 2 )) = Voi(T n SOuo 1 )) = [ 4 A (Ri,R2,d) dr = 5 96793 ,,. 

2 

Proceeding in the same way, one can compute 

Vol(T n S(ui)) = Vol(T n S(u 3 )) = 5.34536 

Vol ((T n S(u^)) \ (SiuJ 1 ) n (5(n 2 ) U 5^ = 2.49982 

Vol ((T n 5(ug x )) \ (5(«3 n (5(1*3) U S(u^))) = 0.70490 • • • 
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Therefore the volume of Q is less than 

36.29366 - 2 • 5.96793 - 2 • 5.34536 - 2.49982 - 2 • 0.70490 = 9.75746 < 9.77029, 
which completes our proof. □ 
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